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Shannon's information entropies in position- and momentum- space and their sum S are calcu- 
lated for various s-p and s-d shell nuclei using a correlated one-body density matrix depending on 
the harmonic oscillator size bo and the short range correlation parameter y which originates from a 
Jastrow correlation function. It is found that the information entropy sum for a nucleus depends 
only on the correlation parameter y through the simple relation S = soa + sia2/ _AsA , where soa, 
sia and XsA depend on the mass number A. A similar approximate expression is also valid for 
the root mean square radius of the nucleus as function of y leading to an approximate expression 
which connects S with the root mean square radius. Finally, we propose a method to determine the 
correlation parameter from the above property of S as well as the linear dependence of S on the 
logarithm of the number of nucleons. 

PACS: 21.10.-k, 89.70,+c 



I. INTRODUCTION 



Information theoretical methods have in recent years played an important role in the study of quantum mechanical 
systems P1J^,|^,[l.|5pP.p|,p|. pi^ , pl]JT^ , [r^ , p^ , p^ ] in two cases: first in the clarification of fundamental concepts of quantum 
mechanics and second in the synthesis of probability densities in position and momentum space. An important step 
was the discovery of an entropic uncertainty relation (!]] which for a three-dimensional system has the form 

S = S r + S k > 3(1 + In tt) - 6.434 (h=l), (1) 

where 

S r = -Jp(r)bxp(r)dr, (2) 

and 

S k = - [n(k)\nn(k)dk (3) 



are the Shannon's information entropies in position- and momentum-space and p(r), n(k) are the density distribution 
(DD) and momentum distribution (MD), respectively, normalized to unity. 

Inequality ([!]) is an information-theoretical uncertainty relation stronger than Heisenberg's and does not depend 
on the unit of length in measuring p(r) and n(k), i.e. the sum S — S r + S k is invariant to uniform scaling of 
coordinates, while the individual entropies S r and S k are not. The physical meaning of S is that it is a measure of 
quantum-mechanical uncertainty and represents the information content of a probability distribution, in our case of 
the nuclear density and momentum distributions. Inequality (|IJ) provides a lower bound for S which is attained for 
Gaussian wave functions. 

Information entropies were employed in the past to study quantum-mechanical systems |l]J|,||,[|,[)|6]j7],||,[| . Recently 
pO| , ^2|Jl4| we studied the position- and momentum-space information entropies S r and S k , respectively, for the densities 
of various systems: the nucleon DD of nuclei, the valence electron DD of atomic clusters and the DD of trapped Bose 
alkali atoms. We found that the same functional form S = a + bhiN for the entropy sum as function of the number of 
particles N holds approximately for the above systems in agreement with Refs. [||]|] for atomic systems. In Ref. 
we found a link of S with the kinetic energy of the system T and a relationship of Shannon's information entropies in 
position-space with an experimental quantity i.e. the root mean square (RMS) radius of nuclei and atomic clusters. In 
Ref. jnj we used another definition of information entropy according to phase-space considerations Q and we derived 
an information-theoretic criterion for the quality of a nuclear DD, i.e. the larger S the better the quality of nuclear 
fj| the DD, the MD and the Shannon's information entropies have been calculated for nuclei using 
uster expansions. The parameters of the various expressions have been determined by least-squares 
fit of the theoretical charge form factor to the experimental one. It was found that the larger the entropy sum the 
smaller the value of x 2 , indicating that the maximal S is a criterion of the quality of a given nuclear model according 



model. In Ref. | 

three different c 
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to the maximum entropy principle. Only two exceptions to that rule were found out of many cases examined. Finally, 
in Ref. []l8f we considered the single particle states of a nucleon in nuclei, a A in hypernuclei and a valence electron 
in atomic clusters. We proposed a connection of S with the energy E of the single particle states. Before proceeding, 
it is appropriate to mention that additional applications of entropy have attracted interest in recent years but 
in a different spirit, in nuclear physics problems, such as in analysis of shell model eigenvectors. However, it is noted 
that in Ref. [[7J the authors define a correlational entropy. This is a von Neumann entropy, which they apply in the 
framework of the nuclear shell model. In our case we use the definition of information entropy according to Shannon 
applied to the density distribution in coordinate, and in momentum space of a nuclear system. 

In the present work we focus on the entropy sum S of a nucleus using the analytical expressions of the DD and 
MD of Refs. The expressions of those distributions have been found for s-p and s-d shell nuclei using 

the factor cluster expansion of Clark and co-workers (22j and Jastrow correlation function which introduces short 
range correlations (SRC). Those expressions depend on the harmonic oscillator (HO) parameter and the correlation 
parameter. We studied the dependence of S on those parameters and we found how this quantity is connected 
with a fundamental quantity, the RMS radius of the nucleus. Finally, we propose a way to determine the correlation 
parameter of the model using the dependence of S on that parameter as well as the linear dependence on the logarithm 
of the number of nucleons. The HO parameter is determined equating the theoretical charge RMS radius of the nucleus 
with the experimental one. 

The paper is organized as follows. In Sec. II, the general definitions related to the correlated DD, MD as well as the 
dependence of the entropy sum on the correlation parameter are given. In Sec. Ill, we study the dependence of the 
RMS radius on the correlation parameter and we give an approximate expression connecting the information entropy 
sum with the RMS radius of the nucleus. In Sec. IV, we present a method for the determination of the correlation 
parameter from the information entropy sum. Finally, in Sec. V, the summary of the present work is given. 



II. CORRELATED ONE-BODY DENSITY OF S-P AND S-D SHELL NUCLEI AND THEIR ENTROPY 



A general expression for the one-body density matrix of N = Z, s-p and s-d shell nuclei was derived in Refs. [£0 21 
using the factor cluster expansion of Ristig, Ter Low and Clark [E2|. That expression depends on the HO parameter 
fro (fro = (^/(muj)) 1 / 2 ), the occupation probabilities of the various states and the correlation parameter y that comes 
from the Jastrow type correlation function 

f(r) = 1 - cxp[-yrf], r b = r/b , (4) 

which introduces short range correlations. The correlation function f(r) for large values of r goes to 1 and goes to 
for r — ► 0. It is obvious that the effect of correlations introduced by the function /(r) becomes large when the 
correlation parameter y becomes small and vice versa. 

The diagonal part of the one-body density matrix is the DD p(r). The Fourier transform of the DD is the form 
factor 



F(q) = J exp[iqr]p(r)dr, (5) 
while the MD n(k) is given by the particular Fourier transform of the one-body density matrix 

n(k) = -^-^ J exp[ik(r - r')]p(r, r')drdr'. (6) 
The expressions of p(r), n(k) and F(q) (in the two body approximation for the cluster expansion) have the forms 

P(r) = jj§tf[0i(r b ) + 2 (r b ,y)], (7) 



n(k) = OtiA-,,) +(/.(/.-,,.//) 



(8) 



and 



F(q) = N d 1 (q b ) + 2 (q b .,y) 



(9) 
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where kb = kbg and = qbo- The terms Oi, 0\ and 0\ come from the one-body term of the cluster expansion of the 

one-body density matrix and the terms O2, O2 and O2 come from the two-body term. Their expressions as well as 
the expression of the normalization factor N a are given in Refs. ]T9| , pO| , pl| . 

From the expressions of p(r) and n(fc) the Shannon's information entropies in position and momentum space and 
their sum S = S r + Sk can be calculated through Eqs. (||) and (0) for p(r) and n(k) normalized to 1. 

For various values of the parameters bo and y and for the N = Z, s-p and s-d shell nuclei: 4 He, 12 C, 16 0, 24 Mg, 
28 Si, 32 S, 36 Ar and 40 Ca we calculated S r , S k and S = S A , treating 24 Mg, 28 Si, 32 S, 36 Ar as Id shell nuclei. 

It is found that, for the above nuclei, S r and Sk depend on both parameters, bo and y, while their sum Sa depends 
only on the correlation parameter y. The calculated values of Sa for the above mentioned nuclei versus 1/y are 
displayed by points in Fig. 1. 

It is seen that Sa is an increasing function both of 1/y and of the number of nucleons A of the nucleus, while Sa 
depends almost linearly on 1/y. For that reason we fitted the numerical values of Sa with the form 

S = S A {y) = s 0A + sia , (10) 

separately for each nucleus, that is the parameters soa, sia and X s a depend on the mass number A of the nucleus. 
The parameter soa is determined from the values of the information entropy sum in the HO case, i.e, 

soa = Sa(oo), 

while the other two parameters are determined by least-squares fit of the values of Sa calculated from Eq. (|l0|) to 
the corresponding ones calculated from Eqs. (|J) and The values of the parameter soa and the best fit values of 
the parameters S\a and \ s a are displayed in Table I, while the values of S = SA(y) found from Eq. (|l^) using the 
above values of the parameters are displayed by lines in Fig. 1. 

It is seen that the simple form of Sa, given by Eq. ([n]), reproduces very well the numerical values of Sa for all 
nuclei considered. Also, there is a systematic trend of the values of the parameters soa, sia and X s a- The parameter 
sqa depends linearly on the logarithm of A. That is expected, as Soa is equal to the information entropy sum in the 
HO case which depends linearly on the logarithm of the number of the nucleons [ ft0| . The parameter X s a has smaller 
values in the closed shell nuclei 4 He, 16 and 40 Ca than in the corresponding neighboring open shell ones. Finally, 
the parameter s±a is almost a monotonic increasing function of A with an exception for the nucleus 40 Ca. 



III. THE DEPENDENCE OF THE INFORMATION ENTROPY ON THE RMS RADIUS 

From the analytical expression of p(r) given by Eq. (^) the analytical expression of the mean square radius of the 
nucleus can be found. It has, in units of the HO size parameter bo, the form 

r\ ee (r 2 )/b 2 = R 1+ R 2 (y). (11) 

The term R\ comes from the one-body term of the one-body density matrix and for the N = Z, s-p and s-d shell 
nuclei has the form 

Rx = j [6ms + 14?72 S + 3077i p + 7077i d ] , (12) 

where r\ n i are the occupation probabilities of the various states. The expression of the two-body term R 2 {y) is given 
in Ref. fisj| . It depends on the parameter y and the occupation probabilities of the various states. 

The calculated values of the RMS radius ri, in units of 60 versus the parameter 1/y for the various s-p and s-d shell 
nuclei are displayed by points in Fig. 2. As in the case of the entropy sum Sa, is an increasing function of 1/y and 
of A. Although we know the analytical expression of we tried to find a simpler expression for it. We found that 
the simple approximate expression 



rb = n(y) = roA + n A ( - ) , (13) 



reproduces well the values of calculated from Eq. (|1 1|) . The parameter tqa is determined from the relation 

roA = r b (oo) = R\ /2 , 
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while the parameters t\a and \ t a are determined by least-squares fit of the values of r b calculated from Eq. ( ]i3| ) to 
the values of r& calculated from Eq. (|TT1). 

The values of the parameter vqa and the best fit values of the parameters ri a a nd X r A for various nuclei are 
displayed in Table I, while the values of r b calculated from the fitting expression ([13|) usingthe above values of the 
parameters are shown by lines in Fig. 2. It is seen that the simple approximate expression (113) reproduces quite well 
the numerical values of r b . As in the case of the parameters soa, s\a and X s a there is a systematic trend in the values 
of the parameters roA, t\a and \ r A- The parameters roA, via and \ t a are increasing functions of A with exceptions 
in the nuclei 36 Ar and 40 Ca for the parameter \ t a- Both of them are considered as Id and 2s closed shell nuclei, 
respectively. 

From the functional forms of Sa(v) and of r b (y) given by Eqs. (fL0|) and (^3[), respectively, the dependence of Sa 
on the RMS radius of the nucleus can be found. That dependence has the form 

SA[r b ) = s 0A + sia (14) 

V r lA J 

Thus, the information entropy sum of a nucleus has been connected with a fundamental quantity, the RMS radius 
of the nucleus measured in units of the HO size b . 



IV. DETERMINATION OF THE CORRELATION PARAMETER FROM THE INFORMATION 

ENTROPY 

In recent works, it has been shown that the information entropy sum of a quantum system (electrons in atoms p|,p| , 
nucleons in nuclei and valence electrons in atomic clusters [ fl2"| and correlated Bose atoms in a harmonic trap Jl4|) 
depends approximately linearly on the logarithm of the number of particles, given by the form 

S = S(A)=a + blnA, (15) 

where a and b depend on the considered system. 

The question that arises is how that property can be used in practice. A possible way is to determine S(A) for two 
nuclei (such that 4 He and 40 Ca) for which there are enough experimental data and then to find a and b of Eq. Jl5| ) 
from the relations 

5(4) 11140-5(40) In 4 

a = , 16 

In40-hi4 ' V ' 



In40-ln4 V ' 

5(4) and 5(40) can be found calculating first the charge form factors of 4 He and 40 Ca from the relation 

Fch(q) = F(q)f cm (q)f p (q)f DF (q), 

where F(q) is the point form factor of the nucleus calculated from Eq. (||), f cm {q) is the Tassie-Barker |23j center-of- 
mass correction and f„(q) and fDF(q) are the correction for the finite proton size and the Darwin- Foldy relativistic 
correction, respectively |24| . 

The parameters bo and y for 4 He and 40 Ca are determined by least-squares fit of the theoretical F c h{q) to the 
experimental, with the constraint that the experimental charge RMS radius is to be reproduced. The values of the 
parameters bo and y as well as the values of x 2 are displayed in Table II. In the same table the values of bo and \ 2 
in the HO case (y = oo) are also shown. In that case bo is determined from the above mentioned constraint. The 
experimental and the theoretical F c h(q), calculated with and without SRC, for the two nuclei are shown in Fig. 3. 

With the values of bo and y determined in the above described way, for the two nuclei 4 He and 40 Ca, we calculated 
the point p(r) and n{k) from Eqs. (^) and (ph and the Shannon's information entropies S r and Sk from Eqs. (0) and 
(||) and their sum, S*(4) and 5(40). Substituting the values of 5(4) and 5(40) into Eqs. ( p^| ) and (jl^), the parameters 
a and b are determined. The calculated values 

a = 5.4029 and b = 0.9360, 

are quite close to the values a = 5.325 and b — 0.858 which have been found in Rcf. JT9J with Skill interaction. 
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Rearranging Eq. ( [To| ) and replacing Sa(v) by S(A) from Eq. (15), we may write 



Using the values of the parameters soa, s\a and X s a given in Table I and the values of S(A) calculated from Eq. 
(|l5|), the correlation parameter y is determined for the other nuclei without any fit to experimental data. The HO 
parameter bo, can be determined now for each nucleus from the relation 



r 



exp ( r ) ^ + r p + 2m 2 c 2 

^ xp is the experimental mean square charge radius of the nucleus and (r 2 ) is the point mean square radius 
calculated from Eq. (Jill). The last three terms of Eq. ( |l9| ) are the corrections due to the spurious center-of-mass 
motion, the finite proton size and the Darwin-Foldy relativistic correction, respectively. 

The values of b and y for the various s-p and s-d shell nuclei determined in the way above described, as well as 
the values of the least-squares errors, in the comparison of the theoretical F c h(q) to the experimental and the entropy 
sum S are displayed in Table II. In the same table the values of bo, \ 2 an d S when SRC are not included (HO case) 
are also shown. From Table II we can see that there is a systematic behavior of the parameter y. The values of y are 
always smaller (that is large correlations) in the closed shell nuclei, 4 He, 16 0, 36 Ar and 40 Ca, than in the neighboring 
open shell ones. It is mentioned that 36 Ar is treated as Id closed shell nucleus. The above behavior indicates that 
there should be a shell effect in the case of closed shell nuclei for the values of the correlation parameter y. Similar 
behavior was found in Ref. p0[ where the determination of the parameters bo and y were made by individual fit of 
the theoretical F c h{q) to the experimental one. From the same table we can see that y is a monotonically increasing 
function of the number of nucleons of the closed shell nuclei. 

The theoretical F c h{q) with and without SRC, as well as the experimental ones for various nuclei have been plotted 
versus the momentum transfer q in Figs. 4 and 5. It is remarkable that without fit to the experimental charge form 
factors, the present method gives good form factors, reproducing the diffraction minima and maxima in the correct 
place. In nearly all cases, the x 2 values found with SRC are better than the corresponding values in the HO case. 
However, the assessment of the quality of the calculated form factors should not be based solely on the values of the 
least-squares errors but also on the fulfillment of the requirement that all the diffraction minima are reproduced in 
the correct place. Thus, comparing the quality of the form factors calculated in the present method with the ones 
calculated with the harmonic oscillator model, we can say that the quality of the form factors are considerably better 
in the former case. All the diffraction minima (even the third one which seems to exist in the experimental data of 
24 Mg, 28 Si and 32 S) are reproduced in the present method while in the harmonic oscillator model they are not. We 
note also that, in the case of the nucleus 36 Ar there are not experimental data for the form factor. The exception 
which appears in 16 O, where the value of x 2 with SRC is worse compared with the value of x 2 without SRC, should 
not be taken as a drawback of the present method. The reason is that there exist many experimental points at low 
momentum transfer where the HO model gives good form factor, while there are a few experimental points at high 
momentum transfers where the present method reproduces these points very well, as well as all the diffraction minima 
and maxima. That can be seen in Fig. 4b. If the experimental points were distributed uniformly, then the values 
of x 2 calculated within the present method would be smaller than the ones calculated within the HO model. Thus, 
we should conclude that even in 16 the theoretical F c h{q) calculated within the present method is better than that 
calculated within the HO model. 

It should be noted that in the above analysis, the open shell nuclei, except 12 C, have been treated as Id shell 
nuclei. However, the same analysis could be made if they were considered as \d-2s shell nuclei, provided that the 
corresponding occupation probabilities of the various states are known. That analysis was made, as an example, for 
36 Ar, assuming that the occupation probabilities of the various states are: rji a = rji p = i]2 S = 1 and rjid = 0.8. The 
values of the parameters bo and y which were found are: bo = 1.7861 fm and y = 7.7658. Thus, the value of the 
parameter y now becomes larger (less correlations) than that of the closed shell nucleus Ca and closer to the values 
of the open shell nuclei 24 Mg and 28 Si (see Table II). 

Finally, comparing the values of the information entropy sum, which were calculated with and without SRC and 
are displayed in Table II, it is seen that the introduction of SRC increases the information entropy sum by 3% to 5% 



bl , h 2 



in agreement with the simple model of SRC used in Ref. 1 1 1 
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V. SUMMARY 



In the present work a systematic study of Shannon's information entropy sum S has been made for various N = Z, 
s-p and s-d shell nuclei using correlated one-body density matrix which depends on the HO size bo and the correlation 
parameter y. 

It is found that, for all the nuclei we have examined, S depends only on y through a simple two-parameter relation. 
A similar approximate expression holds for the RMS radius of the nucleus fb, measured in units of bo- It is found also 
that a simple relation connects S with tv 

From the dependence of S on y and its linear dependence on the logarithm of the number of nucleons of the nucleus, 
the correlation parameter y for a nucleus can be determined, provided that there are enough experimental data for 
two neighboring nuclei. It is mentioned that, usually, the two parameters of the correlated one-body density matrix 
are determined for each nucleus by least-squares fit of the theoretical F c h{q) to the experimental. Within the present 
method, those parameters are determined even in those cases where there are not any experimental data for the 
charge form factor as this has been made for the nucleus 36 Ar. The only experimental data which are used are the 
experimental charge RMS radius of the nucleus, as well as the experimental charge form factors and RMS radii of 
only two nuclei, those of 4 He and 40 Ca. The procedure which is followed to determine the parameters of the model 
is the following: First we determine the parameters a and b of Eq. (15) from the experimental data of 4 He and 40 Ca. 
Thus, the values of S(A) for each nucleus are known. From the value of S(A) for a particular nucleus and from 
Eq. (|l0|) we determine the correlation parameter y corresponding to each nucleus. From that value of y and from 
Eq. (|ll]) the parameter bo is determined so that the theoretical RMS charge radius to be the experimental one. It 
is noted also that, using the sum of the information entropies S r and Sk, the short range correlation parameter has 
been determined indirectly, from the density distribution, as well as from the momentum distribution. This appears 
to be an interesting feature of the present method, since that parameter is usually determined only from the density 
distribution. We would like to mention that "experimental data" for the momentum distribution are not directly 
measured but are obtained by means of y-scaling analysis ]3C| ] and are only known for 4 He and 12 C. 

It would be interesting if the dependence of the information entropy sum on the various parameters could be 
examined with more sophisticated models, such as the density dependent interactions, and if the present method 
could be applied for the determination of a parameter of those models. 
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TABLE I. The values of the parameters s a, «ia and X s a of the information entropy sum Sa of the relation (10), and of 
the parameters roA, tia and X t a of the RMS radius of the relation (13) for various s-p and s-d shell nuclei. 



Nucleus 


S()A 


SlA 




roA 


riA 




\ r A 




4 He 


6.4342 


1.0410 


1.0064 


1.2247 


0.4740 




1.5221 




12 C 


7.5086 


2.1885 


1.1548 


1.4720 


1.5579 




1.8282 




16 o 


7.6069 


2.6464 


1.1529 


1.5 


2.5069 




1.9549 




24 Mg 


8.0933 


3.7445 


1.2390 


1.6330 


4.1480 




2.1115 




28 Si 


8.2096 


4.1641 


1.2548 


1.6691 


5.2681 




2.1959 




32 g 


8.2901 


4.5837 


1.2659 


1.6956 


6.2356 




2.2275 




36 Ar 


8.3490 


4.9578 


1.2681 


1.7159 


6.2862 




2.1758 




40 Ca 


8.4347 


4.7275 


1.2208 


1.7321 


6.5364 




2.0809 




TABLE II. The values of the parameters bo (in 


fm) and y, the x 2 , the RMS charge radius (r 2 ^) 1 ^ 2 


(in fm) and the information 


entropy sum S for various s-p and s- 


d shell nuclei 


. The theoretical RMS charge radii are 


equal to the experimental of Ref. [25]. 


Nucleus 


Case 


bo 


y 


x 2 




(r 2 ch ) 1/2 




S 


4 He 


SRC 


1.2497 


3.7857 


9.40 




1.676 




6.7068 




HO 


1.3335 


oo 


53.46 




1.676 




6.4342 


12 C 


SRC 


1.5617 


7.1294 


153.46 




2.471 




7.7351 




HO 


1.6108 


oo 


181.31 




2.471 




7.5086 


16 o 


SRC 


1.6451 


5.1782 


417.05 




2.730 




8.0044 




HO 


1.7554 


oo 


202.09 




2.730 




7.6069 


24 Mg 


SRC 


1.7609 


7.8711 


221.07 




3.075 




8.3839 




HO 


1.8222 


oo 


226.41 
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FIG. 1. The information entropy sum 
Sa versus the correlation parameter 1/y 
for various s-p and s-d shell nuclei. The 
points correspond to the numerical values 
of Sa(v) and the lines come from the fit- 
ting expression (10). 
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FIG. 2. The RMS radius r b , in units of 
the HO size bo, versus the correlation pa- 
rameter 1/y for various s-p and s-d shell 
nuclei. The points correspond to the val- 
ues of Tb calculating from Eq. (11) and the 
lines come from the fitting expression (13). 
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FIG. 3. The charge form factors of nuclei 4 He (a) and 40 Ca (b). The solid lines correspond to the case when SRC are included 
and the parameters y and bo are determined by least squares fit of the theoretical charge form to the experimental with the 
constraint the calculated RMS charge radius is to be the experimental one. The dot lines correspond to the HO case when bo is 
determined from the experimental RMS charge radius. The experimental points for 4 He are from Ref. [26] and for 40 Ca from 
Ref. [27]. 
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FIG. 4. The charge form factors of nuclei 12 C (a), 16 (b) and 24 Mg (c). The solid lines correspond to the case when SRC 
are included and the parameters y and 60 are determined from Eq. (18) and the experimental RMS charge radius, respectively. 
The dot lines correspond to the HO case when the parameter bo is determined from the experimental RMS charge radius. The 
experimental points for 12 C and 16 are from Ref. [28] and for 24 Mg from Ref. [29]. 




FIG. 5. The charge form factors of nuclei 28 Si (a), 32 S (b) and 36 Ar (c). The various cases are as in Fig. 4. The experimental 
points for 28 Si and 32 S are from Ref. [29]. 
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